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1. Introduction 



In this article, the concentration of the study is on the bivariate extreme distri- 
bution for minima, BEVM, with a specific dependence assumption which lead 
to the development of restricted Pickands' dependence function. We develop 
bivariate extremal models and associated statistical procedures for vector ob- 
servations whose components are subject to an order relationship. The joint 
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structure of the extreme value models for ordered variables are needed to ex- 
tract the information from extreme value data. We restrict our scope of study 
by considering the order of the two variables involved i.e. one variable is greater 
than the second variable. 

We elaborate the theoretical properties of the bivariate extreme value distri- 
bution for minima, BEVM with constraint of the variable X is less than variable 
Y, X < Y. We include the dependence case with X and Y in the form of expo- 
nential margin, but with Frechet margin and the generalized extreme value for 
minima (GEVM) margin also being used in certain places. 

The main difference of this article compared to Nadarajah, (6), is in our 
development of the BEVM model, we are dealing with the restricted logistic 
function and GEVM marginal distribution. We introduce some new terminolo- 
gies and theories related to the modification of dependence functions, A. We 
use Pickands non-parametric estimator of the dependence function to guide us 
in the analysis. 

The structure of this paper is as follows. In Sections 2 and 3 we explain the 
theoretical background of the bivariate extreme value distribution and ordering 
properties of BEVM. Later in Sections 4 and 5, we discuss the asymmetry lo- 
gistic function leading to a new formulation of a restricted logistic dependence 
function and a restricted exponential measure function. We also illustrated some 
features of the dependence function in Section 6. Later in the paper, we also il- 
lustrated how to estimating the boundary constant value c using non-parametric 
estimation of dependence function. At the end of this paper, an application using 
simulation data is also demonstrated. 

2. Bivariate Extreme Value 

Bivariate extreme value theory can also be used to examine the dependence 
between two series of extreme data, X and Y. Let (X, Y) follow the bivariate 
extreme value distribution which arises as the limiting distribution of compo- 
nentwise maxima see (3). He studies the BEV for maxima but as here we are 
interested in minima we will consider the BEVM distribution with exponential 
margins initially. The general bivariate extreme survival distribution function 
for random variables X and Y with exponential margins, is defined as 



G(x, y) — Pr(V > x, Y > y) = cxp[— V(x, y)] = cxp 



-(x + y)A 



( 



x + y 



x 



) 



} 



where the exponential measure function 




with 




(2.1) 
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where the characteristic distribution function H(q), is a finite positive measure 
over [0, 1], and dependence function 

A(u>) = [ max{(l - w)q,w(l - q)}dH(q) 
Jo 

= u (l-q)dH(q) + (l-w) qdH{q) (2.2) 

where (1 — u)q > u>(l — q) if and only if > yz^j, and hence if and only if 
< uj < q < 1. Pickands, (7), defined A(co) as the dependence function of (X, Y). 
From equations (2.1) and (2.2) we have A(0) = A(l) = 1, max(cj, 1 — uS) < 
A(lu) < 1 for < ui < 1 and A(u>) is a convex function in the region < u) < 1. 
The estimation of the Bivariate Extreme Value (BEV) distribution reduces to 
estimation of the A function when we assume that the marginal distributions 
have already been transformed to exponential, see (10) and (9) for more features 
of A(u). 

The joint survivor function G(x, y) is uniquely determined by specification 
of any of V, A and H functions. The stronger dependence between X and Y 
arises when H places mass closer to the centre of [0,1] in which case A(u>) 
is closer to the bounding curve max(cj, 1 — co) for uj G [0,1]. The X and Y 
are totally independent if A(u>) = 1. The X and Y are perfectly dependent 
if A(tu) = max(w, 1 — ui). From equation (2.2) the first derivative of A(ui) is 
A'(uj) = H{uj) — 1, as J Q qdH(q) = 1. We get the relation 

H(w) = A'(ui) + 1 (2.3) 

and clearly h{oS) = H'(u>) = A"(u>). 

Tawn, (10-12) explained in detail about the bivariate extreme distribution 
and more generally the multivariate extreme distribution. Here, we use only the 
logistic model because it is widely used due to it's simple structure. Some details 
about BEV distribution are also available in (1), (2) and (9). 

A widely used non-parametric estimate of A, and here V(x,y) and G(x,y), 
is the Pickands' estimator see (7; 8). Pickands' estimator is defined as follow 

1» = H (2.4) 

with T w , 2 = min . But A p (0) = | + 1 and A p (l) = | j4 1. We use a 

modification of estimator (2.4) proposed by (4; 5) where they suggest replacing 
T u ,i by 



X, Y, 



lux (i-w)y 



(2.5) 



This estimator satisfies A p (0) — | = 1 and A p (l) — i = 1 and A p {u>) > 
max(u, 1 — uS). 
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We want to develop a related statistical methods for vector observations 
whose components are subject to an order restriction, X < Y. We use an asym- 
metric logistic function from (6) and (10) as our dependence function. When 
the order restriction is accounted for we obtain a restricted logistic function, a 
modification from the asymmetric logistic function. We obtained the exponen- 
tial measure function for the restricted case and later we study the features of 
the restricted dependence function. 



3. Ordering Properties of BEVM 



We get the BEVM^, (Tj,, jKj,, (Tj,, ^y) distribution function from the survival 
function of the extreme minima data where 



G(x,y) = P(X > x,Y > y) 
exp • 



1 - U) 



C(x,y) 



— dH( U ) 



with C(x,y) = 1+D \ X:y} , where D{x,y) = ^-Ljjf kr- For x in thc ran S c 
of Y,(x £ Ry), when we take Y = x then for X < Y, 

-dG(x.y) I _ 

Pr(X > x\Y = x) = 9y - . V ~ X = 0, 

where g(x) = — dC ' q£' x ^ , marginal density distribution for Y . As g(x) ^ 0, then 



for Pr(X > x\Y = x) = we need 9G ^ y) \ y=x 

-G(x,y) rC{*,y) 



dG(x,y) 
dy 



as 



1-C(x,y) 



1+1/ tv 

+ 



= 0. However, 

(l-Lu)dH(uj), 



C(x,y) 



When y = x, 

dG(x,y) 



dy 



-G(x,x) 



y=x 



+ 



C(x,x) 



(l-w)dH(ui) 
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So as Pr(X > x\Y = x) = for all x £ R Y holds only when J^ x,x) (l - us)dH(uj) = 0. 
This implies that 

h{ui) = for allw £ [0, C(x, x)\, (3.1) 

for a given a; 6 Ry- As we need this to be zero for all x so h(oj) = for allw £ 
[0, c] where c = max x gfl y C^c, x) or c = j-xg where d = min a . G .fj y . -D(a;, x). 

4. Asymmetric Logistic Function 

We want to develop a bivariate extreme value model and related statistical 
methods for vector observations whose components are subject to an order re- 
striction, X < Y. We use an asymmetric logistic function from (6) and (10) as 
our dependence function. When the order restriction is accounted for wc obtain 
a restricted logistic function, a modification from the asymmetric logistic func- 
tion. Wc obtained the exponential measure function for the restricted case and 
later we study the features of the restricted dependence function. 

Nadarajah, (6), and Tawn, (10), defined the exponential measure function 
for asymmetric logistic model with exponential margins as follows 

V a (x, y) = f max [ujx, (1 - u>)y\ dH a (u>) 
Jo 

with 

h a {oj) = ( s - i)(^ 2 )>(i - w)] s - 2 [(^) s + e s 2 (i - uj) s ^- 2 , (4.i) 

and H a ({0}) =1-0! and H a ({l}) = l-9 2 , with 

i ,i 
u}dH a {uj)dio = (1 — ui)dH a (ui)duj = 1. 
Jo 



But, 



V a (x,y) = {x + y)A a 



, x + v, 

From (10) if h(u>) is from the asymmetric logistic function then we get 
Aa(w) = (1 - 9i )u> + (1 - 6 2 ){l - w) + {(9M S + [0 2 (1 - uj)] s }^ 

A a (oj) = (i - e 1 ) - (i - e 2 ) + [ei^- 1 - e s 2 (i - uy- 1 }^* + e s 2 (i - u,)*]*- 1 . 

From equation (2.3) where H{ui) = A'(lo) + 1 we get 

H a (uj) = 1 - 6 1 + 9 2 + [^cj 5 - 1 - 05(1 - w)"][0;o;' + e 8 2 (l - w) 8 ]*" 1 , <oj<1 
and then 

H a (0+) = (2 - fr) - (1 - a ) + (-0|) (O^- 1 = 1-6,, 

and 

fr a (i_) = 2 - 0! - (i - 2 ) + e{e\- s = i + 2 . 
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5. Getting a Modification Dependence Function 

We use the notation N for the logistic function used by (6) which is defined on 
(a, b) , with 0<a<^,^<6<f, and RL for our restricted function which is 
defined on (c, 1), with < c < \. 
Using Theorem 5.5 from (6) gives 

M«) = ( s - !)( & - a){a(3) s [(w - c)(6 - w)] s " 2 

x [a s (w-a) s +/3 s (6-u;) s ] = ~ 2 , we (a,6),s>l (5.1) 

as a valid density of the measure Hn on [a, b] with _ffjv({a}) = 7i and iijv({fr}) = 
72 and satisfying equation (2.3). Given a < i, 6 > i, 71,72 > are such that 
71 + 72 < 2, 72 < and 71 < T=JT. and a = 2fe - 1 and /? = 1 - 2a. 

From equation (3.1), /i(w) = for w e [0,c] then the solution for is 

Vjv(a;, y) = [ max [wa;, (1 — hj\r(u))dcj = J max [wi, (1 — }in(uj)cLj. 

(5.2) 







Using equation (5.1) with a = c, b = 1, a = 1, /3 = 1 — 2c, w e (c, 1) we focus 
on the restricted function so 

= (s - 1)(1 - c)(l - 2c) s [(w - c)(l - w)] s " 2 

x [(w - c) s + (1 - 2c) s (l - cj) s ]i- 2 . (5.3) 

Comparing equations (4.1), (5.1) and (5.3), we get 6\ = a = 1 and 6> 2 = /? = 1 — 2c 
with 71 = 72 = 0. 

We transform the function with w = y5f => w = c + (1 — c)tt), < w < 1, 
and J^j = 1 — c, we get then the exponential measure for the restricted logistic 
function as follows 



f 1 h, 
V RL (x,y)= J max{[c+ (1 - c)w]x, [(1 - c)(l - w)]y} -y 



h a {w) 

■aw 



max 



11 



1 - c 



+ w 1 x, (1 — w)y 



h a (w)dw 



with h a (w) defined in equation (4.1), with 6% = 1 and #2 = 1 — 2c. It is easy to 
show that 



V RL {x,y) 



for 0<^<c 



({[(l-c)y-ca;] s + (l-2c)^}* + cx) for c < ^ < 1. 

(5.4) 



The proof is provided at the end of this paper. 

Now we identify the characteristics of new dependence function, A m, , of the 
restricted logistic dependence function by plotting Arl(w) versus w. 
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0.0 0.2 0.4 0.6 0.8 1.0 



Fig 1. We fix c = 0.25 but with varieties of s = 1 (solid line), 1.5 (dashed line) and 2.5 
(dotted line) values. 




Fig 2. We fix s = 1 but with varieties of c = 0.1 (solid line), c = 0.3 (dashed line) and 
c = 0.45 (dotted line) values. 

In Figure 1 wc fixed c = 0.25 but with three values of s, where s = 1.00, 1.50 
and 2.50. When X < Y with c = 0.25, the X and Y are as weakly dependent 
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as possible if 

. . , f l-w for < w < 0.25 
A RL {w) = ^ , w+ 2 for .25 < W <1. 

By increasing the values of ,s to 1.50 and 2.50 we get a stronger dependence 
between X and Y in which A(w) is closer to the bounding curve max(w, 1 — w) 
for w e (0.25,1]. 

In Figure 2 we study the features when s = 1, i.e. the weakest possible 
dependence between X and Y given the value of c. By increasing the values of 
c, -A(iu) values are given in Figure 2. 



6. Extension Cases for Ordering 

In this section, we will derive the dependence function Am for the Nadarajah 
model when hpf(uj) = 0,uj e (0, c) then we propose the dependence functions, 
An for when hjsr is only non-zero in [0, c] then in [ci, C2]. 

Proposition 6.1. Let An(u>) be the dependence function of (X,Y) when the 
measure h^ is Nadarajah logistic function with h^{u>) = for to G (0, c), < 
c < i. then An{oj) function is 

fl-w J /0<w<c, 

j£(l^) + _l_ [(c ,_ c)S + (1 _ 2c)S(1 _ w)1 i i/c<w <l. 

1 < s. 

Proof 

For < u> < c, we get 

A' N (w) = -1, A^(w) = = /uv(o;). 
For c < a; < 1, we get as follow 

(s - 1)(1 - c)(l - 2c) s [(w - c)(l - uj)] s - 2 [(w - c) s + (1 - 2c) s (l - . 

Compared with the density measure of equation (5.3), then ^4^(0;) = h^(uj). 

When X < Y, h(uj) = for < cj < c, where < c < 5. But this is not the 
only stochastic ordering of interest in general, other ordcrings might be studied. 
Possible ordering as X < Y* where Y* can be any function such as 1/Y, logF, 
Y 2 etc. We might also be interested in X > Y. We propose a possible result for 
new dependence functions for any ordering of X < Y* where the dependence is 
of a logistic type such that the following two cases follow 

1. If h^iio) = for c < u> < 1, where i < c < 1. 

2. If /ijv(w) = for < w < c\ and C2 < w < 1 where < ci < | and 

i < C 2 < 1. 
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For any possible relationship between X and Y, that leads to ft.jv(w) = for 
c < uj < 1, where | < c < 1 we get 

Vn(x, y) = I max (wi, (1 — ui)y) h^{uj)duj. 
Jo 

Proposition 6.2. Let An(oj) be a dependence function of (X, Y) when the 
measure h^ is Nadarajah logistic dependence function with hjsr(u>) = for c < 
uj < 1, with | < c < 1 then An(w) function can be written as 



A n (lu) 



' Q-^u + 1 p c _ ^s^s + (c _ ifO<OJ<C, 

uj if c < uj < 1, 



with 1 < s. 



When /liv(w) = for < u> < c\ and c 2 < w < 1 where < ci < i and 
5 < c 2 < 1. 

Vjsr(x,y) = / max (wx, (1 — oj)y) hN(u>)du>. 

Proposition 6.3. Let An(oj) be a dependence function of (X, Y) when the 
measure is Nadarajah logistic dependence function with h^{uj) = for < 
to < c\ and c 2 < uj < 1 with 0<Ci<|<C2<l, then An(uj) function is 

!1 — uj if < uj < a , 

(1 T 2 -r ci) + £ fer 1 + w/ H *M 
UJ if C-2 < UJ < 1 , 

where W{uj) = [(2c 2 - 1) s (uj - a) s + (1 - 2ci) s (c 2 - w) s ] • with 1 < s. 

The proof of the latter two propositions similiar to Proposition 6.1. We show 
the dependence function by using the Proposition 6.3. In Figure 3 we fixed 
Ci = 0.25 and c 2 = 0.75 but with three values of s, where s = 1, 2 and 5. When 
s = 1 , we get the weakest possible dependence for given c\ , c 2 when 

!1 — uj if < uj < a , 
0.75 ifci<w<c 2 , 
w if c 2 < cj < 1. 

By increasing the values of s from 1 to 2 and 5 we get a stronger dependence 
between X and Y in which An(uj) is closer to the bounding curve max(w, 1 — a;) 
for uj e (0.25, 0.75). In Figure 4 we show A n (uj) for fixed ci = 0.10 and c 2 = 0.60, 
s = 1, 2 and 5. 
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~ i i i 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 



Fig 3. We fix a = 0.25 and b = 0.75 but with varieties of s = 1 (solid line), s = 2 (dashed 
line) and s = 5 (dotted line) values. 




~ i 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 



Fig 4. We fix a = 0.10 and b = 0.60 but with varieties of s = 1 (solid line), s = 2 (dashed 
line) and s = 5 (dotted line) values. 

7. Estimating c Using a Non-parametric Estimator of Dependence 
Function 

In this section we show how to estimate c using Pickands' non-parametric es- 
timator of the dependence function. We show also that the estimated c value 
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obtained from Pickands' estimator is always greater than the true value of c. As 
X < Y we know that 

A(w) = l-w for < w < c (7.1) 

for some < c < 0.5. We wish to estimate c using Pickands' estimator defined 
in equation (2.5), the modified Pickands estimator by (5) 

A p (w) 



IXi min (i?'7r^) 



For AJw) = 1 — w for given < w < 0.5 then -4i > - — Vi, i = 1, . . . , n 

yx ' wX (l—w)Y 

Let Wi = x X _i Y - 1 tncn 

^ ^ Yi ^ Wi ^ l-Wj w . . 1 

> = =>• > =, vi, i = 1, ... ,11. 

wX (l-w)Y wX (l-w)Y 

=>■ Wi > 



(i - + wx 

= = Vi, i = 1, . . . , n. 

As n — > co X — > 1 and V — > 1 so — =4 r — ► This leads to our conclusion 

i+=(^-i) 

that 

^ Wi > w Vi, i = 1, . . . , n. 
So if A(w) satisfies equation (7.1), modified estimator satisfies 

c < ^Pickands = min(Wi. ■ • • , Wn)- 
We recognise this is biased as c < cpickands always. 



8. Simulation study 

We simulate 50 repetitions from (Z x ,Z y ) ~ BEVM(s, fj, x (t),a x , n y (t),a y ,£) 
where s = 2, /i x (t) = 100 - 40t, ^(t) = 150 - 40i, cr^ = 4, cr y = 2 and £ = 0.2. 
This distribution satisfies Condition 1 of Theorem 2 in Chapter 6 with the 
marginal distribution for Z x is QEVM.(n x (t),o x ,£) and Z y is GEVM(/^ y (£), cr y ,£). 
These data are non-stationary with a trend in mean i.e. fi x (t) = 100 — 40t and 
H y (t) = 150 — 40i which later are estimated using penalized log-likelihood, i.e. 
g x and g y respectively. 

Figure 5 (left) shows realisations of Z x and Z y against time showing clearly 
that Z x < Z y for all realisations. Figure 5 (right) shows that there is a depen- 
dence relationship for Z x and Z y after de-trending (removing the trend to give 
IID residuals). 
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Fig 5. Left plot is a scatter-plot for marginal data of Z x (black circle) with trend fi x (t) = 100 — 40t (black line) and g x (black dashed line); and Z y 
(red circles) with trend [iy(t) = 150 — 40i (red line) and g y (red dashed line). Right plot is Z x — g x versus Z y — g y . 
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Table 1 

The convergence values for parameters s, n x , fi y and £ from the simulation data. 



s 


0~x 


0~y 


t 


1.68 


3.71 


1.85 


0.214 




200 400 600 800 1000 
iterations 



200 400 600 800 1000 
iterations 



Fig 6. The plot for converged parameter values for a x , cr y , £ and s using penalised log- 
likelihood with smoothing parameter X x = 1000 and X y = 1000 for 1000 iterations. 



From Figure 6 we see the convergence for iterative estimates of the parameters 
s, (T x , o x and £, with final estimates tabulated in Table 1. We have the g x and 
g y represented by dashed lines in Figure 5 (left). 

Figure 7 shows that there is strong dependence between Z x and Z y as all 
dependence functions are closer to the bounding curve max(w, 1 — w) for w € 
[c, 1] with values for each c for dependence functions. 

From the working of the theory, the density measure fiRL is defined on 
c = 0.030 < Cpickands ~ 0.078 for true marginal parameter values and c — 0.052 < 
Cpickands = 0.070 for estimated marginal parameter values from the simulation 
data. 

Before assessing the model fit using P-P and Q-Q plots we need to transform 
the data to Exponential margins using equation (7.37) discussed in Subsec- 
tion 7.5.2. By removing the marginal aspect in the data, the P-P and Q-Q plots 
in Figure 8 show that the data give a good fit to the BEVM of the restricted 
logistic function for the simulation data in Figure 5 when there is assumed to 
be dependence. 
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Fig 7. Simulation data: Pickands dependence function estimates, A(w) based on annual min- 
ima of Figure 5. N on- parametric estimates by Pickands using true parameter values ( red line ) 
and using estimated parameter values (dashed red line). Parametric estimates for restricted 
logistic function using true parameter values (black line ) and using estimated parameter val- 
ues (dashed black line). 



9. Discussion 



When X and Y are ordered with the special constraint of stochastic ordering 
X < Y, it imposes som climitations to the model of ordered bivariate extreme 
values. The logistics function defined by (6, Nadarajah) and (10, Tawn) was 
used to obtain a modify function called a restricted logistic function to model 
the dependence between variables. We proposed also some extensions of the 
dependence function of Nadarajah more generally. One of the important results 
is that when ordering exist, i.e. X < Y, we found that the density measure 
h(w) = 0, for w £ [0, c = maxjegfly C(x, x)}. This led to some other result, i.e. 
the restricted logistic dependence function (a modification version of asymetrical 
logistic function), restricted exponential function, Vrl(x,j/). 

We showed through a simulation study that the ordered bivariate moel fits 
well when the constraint X < Y is imposed. The simulation study with selection 
of parallel linear trends for X and Y showed that this approach can be applied 
to real data. 
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Fig 8. The P-P (left) and Q-Q (right) plots for the simulation data of(Z x ,Z y ) ~ BEVM(s = 
2,fj, x (t) = 100 - 40t,o x = 4,/ij,(t) = 150 - 40t, <7 H = 2,£ = 0.2) w/iic/i followed to Condition 
1 of Theorem 2 in Chapter 6 with 50 repetitions. 



Proof of equation (5.4). 

The measure distribution function, H a (w), is derived as follows 



H a (w)= / h a {t)dt 
Jo 

= 1 - 2c + [w 3 ' 1 - (1 - 2c) s (l - wy- 1 }^ + (1 - 2c) s (l - w)^- 1 . 
We get the integral of H a (t) for < t < w and < t < 1 as follows 

H a (t)dt = -(1 - 2c) (1 - w) + [w s + (1 - 2c) s (l - w) s ] • , 



/ H a (t)dt=l, 
Jo 



and £T a (0) = and #„(1) = 2(1 - c). 

The exponential measure function for restricted logistic function we simplify 
as follows 

Vrl{x,v) =y [ (1 - w)h a (w)dw + / h a {w)dw + x f wh a {w)dw, 

Jo l — c 7 D+ 

where D + = max(D,0) with < D + < 1 and 

x + j/ \1 — cj x + y \1 — cj x + y 1 — c \x + y J 



imsart-ejs ver. 2008/08/29 file: ejs_2009_358.tex date: January 12, 2009 



Mohd Bakri Adam/ Bivariate Extreme Modelling 



16 



When c < —j— < 1 then D + = D, the derivation of Vrl is using integration by 
parts, letting u = 1 — w and dv = h a (w)dw du — ~dw and v = H a (w), we 
get 

/ (l-w)h a (w)dw = L(D + ) + M(D + )P{D + )(l-D+) (9.1) 
Jo 

where 

L(D+) = [D\ + (1 - 2c) s (l - £+) s ] i , (9.2) 

M(Z>+) = + (1 - 2c)'(l - £+) s ] , (9.3) 
P(D+) = D^ 1 - (1 - 2c) s (l - D+y- 1 . (9.4) 

Letting u — w and cfo = h a (w)dw du = dw and i> = H a (w), 



[ wh a (w)dw = L{D+) - D + M{D+)P{D+) 
Jd + 

f h a {w)dw = 1 - M(D+)P(D + ). 



(9.5) 
(9.6) 



Using equations (9.1), (9.2), (9.4), (9.3), (9.5) and (9.6) then 

V RL (x,y)=(y + x)L(D + ) + - -, 

(1-c) 

as 

cx 

yM(D + )P(D + ) - - -M(D + )P(D + ) = yD + M(D + )P(D + ) + xD + M(D + )P(D, 

(1-c) 

When y = then D + = 0. 

Then the exponential measure function is 

Vrl(x,0) = x. 

When x = then D + = 1. Then the exponential measure function is Vrl(0, y) = 

y- 

i 

< c then D + = we get 



When c = 0, then L> + = ^ we get V RL {x,y) = {y s + x s ) s . When < 



x+y 

V RL (x,y) = 2cx+(l-2c)x = x. (9.7) 

Checking for c = then when x = 0, we get V(0, y) = y. When y = 0, we 
get V(x, 0) = x. Another way of writing the exponential measure function is by 
including the dependence function, A, 

V RL (x, y) = (x + y) A RL ( ) . (9.8) 
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From equation (9.7), for < w — -H— < c 



Arl (w) 



x+y 

x = (x + y) A RL 

X 



x + y 



x + y 



= 1 — w. 



When c < w = < 1 => D + = D from equation (9.2) and relation with 
equation (9.8) we obtain 



A 



RL 



H = — !— {c(l - to) + l(w - c) s + (1 - 2c) s (l - w) s }^} 



Checking when w = c then 

Arl (c) = 7^— [c(l - c) + (1 - 2c)(l - c)] = c + 1 - 2c = 1 - c. 

1 — c 

Considering all possible values for D + the generalized form of the exponential 
measure for restricted logistic is 



V RL {x,v) 



i=j ({[(1 - c)y - cx] s + (1 - 2c) s a; s }* + ex) 



for < < c 

for c < < 1. 

x+y 



Using equation (9.8) we can easily get the restricted dependence function of 
Arl{x, y) where for < w < c 

Arl(x,v) = 1-oj 

and for c < u> < 1, we get 

A RL {x,y) = _i_ | c (l - a;) + [(w-c)* + (l-2c)° + (1 - 2c)« (1 - w)-]"} . 

From (3) the relationship between an exponential measure, V, and a density 
measure function, h, involving two variables of Frechet margins denoted by "F" 
where x — — and y = — is 



2 Vrl {x F ,yF) 
dxpdyF 



(x F + y F f 



IRL 



%F 



XF + Uf 



By changing the equation (5.4) to Frechet margins we get Vrl(xf, Vf) is equal 
to 



i 

XF 



for < — 2£— < c 



RL 



1 

1-c 



1 — c c 

VF XF 



for c < 



XF+y f 
(9.9) 



< 1. 
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0. As y F -> oo 



h RL (0) = 0. We get for 0< x 
d 2 V RL (x F ,yF) 



xf+Vf 



< C 



1 



dxpdyF 



(x F + yF) 3 



IRL 



VF 



x F + yF 



= 0. 



From equation (9.9) checking for c < ^ - < 1, we get 



d 2 VRL(x F ,yF) 
dxpdyF 



1 



{y F x F ) 3 
1 - c 



(l-c)(*-l)(l-2c)« 



1 / 1 - c 
J'f 



2/F 



-— +(l-2c)» — 



1 



Xf 



Vf 

-2 



IS 



we know 



c 



s-2 



(9.10) 



= - (s - 1)(1 - c)(l - 2c) s {y F (z F (l - c) - q/ F )r 
x [(x F (l - c) - cj/ F ) s + (1 - 2c) s y F ]^ 2 . 



■V F 



X F + yF 



( S -l)(l-c)(l-2c) s 



XF + DF 



VF 



XF + VF 



xf 



1 



•IF + VF 

3 



-c +(l-2c) s 



2/F 



x F + yp 



XF + VF 



(s - 1)(1 - c)(l - 2c) s {j/f [a*(l - c) - cy F ]} 
x {[z F (l - c) - C2/F ] s + (1 - 2c) s z^}^ 2 
(xf + Vf) 



s-2 



^d 2 V RL {x F ,yF) 



dxdy 



d 2 Vrl {xf i ?/f ) 
dxpdyF 



1 



(a?F + Vf) 3 '" RL \xf + Vf 



xf 
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